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For a single membrane of stiffness k fluctuating between 
two planar walls of distance d, we calculate analytically the 
pressure law 
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The prefactor tt 2 /128 ~ 0.077115 . . . is in very good agree- 
ment with results from Monte Carlo simulations 0.079±0.002. 

1. A stack of n parallel, thermally fluctuating mem- 
branes exerts upon the enclosing planar walls a pressure 
which depends on the stiffness K and the temperature T 
as follows 
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where ks is Boltzmann's constant and d the distance 
between the walls (see Fig. 1). 
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FIG. 1. Membrane fluctuating between walls of distance d, 
exerting a pressure p. 

This law, first deduced from dimensional considera- 
tions by Helfrich Q , is of fundamental importance in the 
statistical mechanics of membranes just as the ideal gas 
law pV — N ksT in the statistical mechanics of point 
particles. We would therefore like to know the size of the 
prefactor, the stack constant a n as accurately as possi- 
ble. So far, its value was determined only by extensive 
Monte Carlo simulations as being 
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For a single membrane, the following value was found 



ai = 0.079 ± 0.002. 



(4) 



So far, there exists no analytic theory to explain these 
values. 

The purpose of this note is to fill this gap for the con- 
stant ai, by calculating analytically the pressure of a 
single membrane between parallel walls. The theoretical 
tool for this has only recently become available: A strong- 
coupling theory developed originally in quantum mechan- 
ics pL was extended successfully to quantum field theo- 
ries [pi, where it has been used to obtain extremely accu- 
rate values for the critical exponents of 0(n)-symmetric 
scalar fields with (/^-interactions 0. 

2. Strong-coupling theory gives direct access to the 
large-g behavior of divergent truncated power series ex- 
pansions of the type 
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The g — > oo -limit of Jn (g), to be denoted by /jy, is ob- 
tained by setting Q = eg 1 / 9 and optimizing the function 



fN(c)=g 1/ "fN(c)=g 
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where 



b N k = 
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is the binomial expansion of (1 — l)^ 1 kq ^ 2 truncated 
after the (N — fc)th term. Optimizing means extremizing 
/jv(c) in c or, if an extremum does not exist, extremizing 
the derivative f' N (c). 

3. We apply this theory to a membrane between walls 
by proceeding as follows. The partition function of the 
membrane is given by the functional integral 

Z = J p u(2 ;)exp J d 2 x[dMx)?} ^e- A f/ k ° T , 

(8) 

where u[x) is a vertical displacement field of the mem- 
brane fluctuating between horizontal walls at u = —d/2 
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and d/2. The quantities A and / are the wall area and the 
free energy per unit area, respectively. Such a restriction 
of a field is hard to treat analytically. 

We therefore perform a transformation which maps the 
interval u S (—d/2, d/2) to an infinite </?-axis, 
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and add to the fluctuation energy E in the exponent of 
(^J) a potential energy which keeps the membrane be- 
tween — d/2 and d/2 (Poschl- Teller potential): 



E pot = E pot + E int = 



d 2 x m A 4> 2 (u(x)) 
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d 2 x < u 2 (x) +^£k 
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with expansion coefficients £2, £3, £4, . 
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FIG. 2. Smooth Potential replacing box walls 

The potential energy per area is plotted in Fig. 1. Its 
presence destroys the simple scaling properties of the 
partition function (^J), which depends only on the di- 
mensionless variable nd 2 /ksT . The new partition func- 
tion Z associated with the modified energy E + E pot has 
an additional dependence on the dimensionless variable 
g = ir 2 /m 2 d 2 . The original hard- wall system is obtained 
in the strong-coupling limit g — > 00. 

In the opposite limit where g goes to zero, the energy 
E + E pot becomes harmonic, 

E o = ^J d 2 x{[d 2 u(x)] 2 + m 4 u 2 (x)}, (12) 

leading to a partition function 

Z = e -i Trl °s( a4+m4 ) = const x e~% m \ (13) 

where A is the area of the walls. 

For a finite distance d, the interaction energy E lnt is 
treated perturbatively order by order in g, expanding the 
exponential e~ E l k sT m a p 0wer series, and each power 
in a sum of all pair contractions. These are pictured 



by loop diagrams whose lines represent the correlation 
function 
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The free energy density / = — ksTA~ x \ogZ is obtained 
from all connected loop diagrams. For simplicity, we shall 
use natural units with n/ksT = 1. 

The lowest contribution to the free energy density 
comes from the expectation value of the u 4 -interaction 



or the loop diagram 3 CO which is of the order l/d 2 
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Together with the exponent in (|13j ) , we thus obtain first- 
order free energy density 
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Continuing the perturbation expansion, yields an expan- 
sion of the general form 
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where 02, ... , aN are dimensionless numbers. By compar- 



ison with (tJ) wc identify p ■■ 
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The function /jv(c) of Eq. (|5|) describing the limiting 
large-g behavior is obtained by setting f2 = c7r 2 /2o? 2 , and 
reads 



a N , N 
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According to the above-described strong-coupling theory, 
we must optimize the expression Jn{c) in parentheses. 
Since the second term does not contain c, we separate 
this term out, and write 



Mc)^^(c)^+(^ + ^- 
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with only the remainder A/at(c) to be optimized. Let 
A/iv be ist optimal value. If we know only 02, we find 
the approximation A/| = a/3c12/16. Ignoring Af N for a 
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moment, the first term in yields the lowest estimate 
for the free energy density of the original system 



fx = 

implying a pressure law 
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By comparison with the general pressure law 
tify the prefactor as being 
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Without the prefactor factor 1/2, this would agree per- 
fectly with the Monte Carlo value (Q). Thus we expect 
the contribution of Af^ for N —> oo to be equal or al- 
most equal to 1/64. 

The calculation of the higher-order terms 02, 03, ... 
is tedious, and will be presented in a separate detailed 
publication |7j] . In this note we shall circumvent it by ex- 
ploiting a close relationship of the present problem with 
a closely analogous exactly solvable one, which may be 
treated in precisely the same way: The euclidean ver- 
sion of a quantum-mechanical point particle in a one- 
dimensional box u G (—d/2,d/2). 

4. The partition function of a particle in a box is 



Z = / Due 



-( K /2k B T) f dx(du) 2 = -Af/kt 
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The quantum-mechanical ground state energy of this sys- 
tem is exactly known: {ksTj k)tt 2 /2d 2 , corresponding to 
a free energy density 
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The path integral ( p4| ) may now be treated as before, i.e., 
we transform u to (f via IM), and separate the field energy 
into a Gaussian energy (in natural units) 



E "='j I dx{[d 2 u{x)] 2 + m 4 u 2 {x)} 



(26) 



and an interaction energy which looks the same as ( |10[ ), 
except that the integration J d 2 x runs now only over one 
dimension, J dx. 

The first-order contribution to the free energy density 
is now (in natural units with k/UbT = 1) 
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with the pair expectation 
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(28) 



leading to a first-order free energy density 
m 2 j_^2 

h ~ ~ + 2d 2 "' 
and a full perturbation expansion of the form 
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From this we find the function /jv(c) defined in Eq. (||) 
governing the strong-coupling limit d — > by setting Q — 



in 
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Here the first term yields the lowest approximation 
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which is precisely half the exact result. Thus we con- 
clude that the optimal value of the neglected expression 
Af;v(c) must be once more equal to 1/4 in the limit 
N —> 00. In order to see how this happens, we extend the 
Bender- Wu recursion relation for the_perturbation coef- 
ficients of the anharmonic oscillator ||. It yields for the 
ground state energy an expansion 
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Inserting the coefficients (|ll]) we find 0,2, a&, a§ . . . : 
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16' 256' 2048' 65536' 524288' 8388608'"'' ^ ' 
whereas the odd coefficinets 03,05,07,... vanish. To 
account for this fact, we resum the series containing only 
the even terms 

c 
4 

taking the coefficients bf of Eq. ^ with the parameters 
p — 1, q = 2. This yields the functions A/jv(c) plotted 
in Fig. 3 and listed in Table 1. 
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FIG. 3. Plots of the functions Af N (c) of Eq. (|3J), all being 
optimal exactly at c = 1 with A/^=l/4. 

For all A/jy(c), optimization yields a strong-coupling 
value A/jy equal to 1/4, thus raising the initial value 1/4 
in (32) to the correct final value 1/2. 

5. To exploit this property of a particle in a box for the 
system at hand, the membrane between walls, we make 
the following crucial observation: The Feynman integrals 
determining the first two terms in the free energy den- 
sities in Eq. ( p"7j ) for a membrane and in Eq. (^9|) for 
a particle are related to each other by a simple trans- 
formation of the integration variables. The membrane 
integrals 
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go over into those of the particle in the box 



J ^ log(q 2 + m 4 ) = 
by the transformation 
k 2 ->o. 
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Thus, if we multiply each loop integral by a factor 1/4, 
we find immediately the free energy density f± of the 
membrane in Eq. (|l7]) from that of the particle in the 
box in Eq. @. 

But the analogy carries further: By differentiation fl36| ) 
and ([I?]) with respect to to 2 , we see that also all Feyn- 
man integrals J d 2 k / (2ir) 2 (k 4 + to 4 )" are related to the 
J dq/(2Tt){q 2 + m 4 ) v by the same factor 1/4. This prop- 
erty has the consequence that most of the connected loop 
diagrams contributing to the perturbation expansion of 
the free energy density, shown in Fig. ^]up to five loops, 
are related by a factor (1/4) L , where L is the number 
of loops. In particular, all such diagrams coincide which 
are ususally summed in the Hartree-Fock approximation 
(chain diagrams, daisy diagrams, etc.). 
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Only the topologicall more involved diagrams 3-1, 4- 
1, 4-2,4-5, 5-2,5-3,5-5,5-6,5-7,5.11,5-12,5-15 in Fig. | 
do not follow this pattern. For a particle in a box, we can 
easily calculate the associated Feynman integrals in x- 
space as described in Chapter 3 of Ref . || , and find that 
they contribute less than 5% to the sum of all diagrams 
at each loop level. This implies that the correspond- 
ing results for the membrane between walls will differ at 
most by this relative amount from those for the particle 
in the box. We therefore conclude that since the optimal 
value of A/jv(c) in Eq. ([32|) doubles the initial value for 
N — » 00, the analogous function for the membrane be- 
tween walls in Eq. (20) will double approximately. For 
the quantitative deviations see the forthcoming publica- 
tion |0| . A precise doubling of the result leads to a 
very good agreement with the Monte Carlo number (g). 

6. The alert reader will have noted that the field trans- 
formation (^) is rather special. We may, for instance, 
chose any mapping 
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(38) 



which has a doubled coefficient of ip 3 with respect to the 
expansio n (H ). As a consequence, the functions /jv(c) in 
( p0| ) and (p2) would have a doubled first term. Since this 
would be the correct final value, the remaining functions 
A/jy(c) would have to converge to a vanishing optimal 
value for N — ► 00 (in the particle case exactly, in the 
membrane case approximately). To reach this goal, the 
coefficients W4, wq, . . . in ( p8[ ) can be chosen rather ar- 
bitrarily, although there are a few convenient ways for 
which the speed of convergence is fast. A preferred choice 
is one in which all coefficients 0,2, 03, 03, . . . of the per- 
turbation expansion vanishes for a particle in a box. This 
and other possibilities will be studied separately [PUTT 
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TABLE I. The functions A/jv(c) of Eq. @, and their 
optimal values A/^. 
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FIG. 4. Vacuum diagrams up to five loops. 
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